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ABSTRACT 


The present paper is concerned with using Bayesian and E-Bayesian method of estimation to find estimates for the shape 
parameter of Exponentiated Inverted Weibull distribution. These estimators are derived by using different loss functions. In this 
paper, Bayesian estimates are derived by using informative prior. 
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1. Introduction 

The Inveted Weibull distribution is one of the most widely used distribution for analyzing life time data. It is found to be useful in 
diverse fields ranging from Engineering to Medical Science. The two parameter Exponentiated Inverted Weibull distribution 
(EIWD) has been proposed by Flaih et al (2012) and is a generalization to the Inverted Weibull distribution through adding a new 
shape parameter 8 € R* by exponentiation to distribution function F. The cumulative distribution function (c. d. f ) of EIWD is 
given by 


Fo (x) = (e**)" ; x, 6,8 >0 


which is simply the Oth power of the distribution function of the Standard Inverted Weibull distribution. Here 8 and B are the 
shape parameters of EIWD. 


The probability density function (p.d.f) of EIWD with two shape parameters 0 and f is given by 
f((x,0,B) = OB x-PtDe®xX? .y 5 0,0>0,8>0 (5.1.1) 


The Exponentiated Inverted Weibull distribution also has a physical interpretation. If there are m-components in a parallel system 
and the life times of the components are independent and identically distributed (i.i.d) as exponentiated inverted weibull distribution, 
then the system lifetime variable also has exponentiated inverted weibull distribution. Many characteristics of Exponentiated 
Weibull Distribution has been found in the literature. Soland (1968) discusses the Bayesian analysis of the Weibull Process. 
Mudholka et al (1995) introduced the exponentiated weibull distribution as a generalization of the standard weibull distribution and 
as a suitable model for bus motor failure time data. Mudholkar and Huston (1996) applied the exponentiated weibull distribution to 
the flood data with some properties. Ahmad et al (2015) discusses the Bayesian estimator of shape parameter of exponentiated 
inverted weibull distribution under different loss function namely square error loss function, entropy loss function and precautionary 
loss function. El-Din et al (2014) worked on statistical inference and prediction for the inverse weibull distribution based on record 
data. Aljouharah (2013) estimates the parameters of ET(WD under type-II censoring. 


Bayesian estimation approach has received great attention by the researchers and they had proved that Bayes estimate perform 
better than classical estimators. A new approach of Bayesian estimation called Expected Bayesian or E-Bayesian method has 
been introduced by Han (2007). He obtained the E- Bayes estimate of failure probability by considering quadratic loss function 
and discussed the properties of E-Bayes estimate and showed that E-Bayes estimate is efficient and easy to operate. [Okasha (2012)] 
constructed Bayesian and E- Bayesian method for estimating the scale parameter, reliability and hazard function of the Weibull 
Distribution under type-II censored samples by considering the squared error loss function. [Jaheen and Okasha (2011)] compared 
the Bayesian and E- Bayesian estimator for the parameters and reliability function of the Burr — XII model under type-II censoring 
by considering the squared error and linex loss function. They have shown that the performance of E-Bayes estimate is better. 
[Azimi et al, (2013)] studied E-Bayesian estimation of generalized half logistic progressive type-ii censored data. [Reyad et al, 
(2015)] develop E-Bayesian estimators of parameter of the Gumbel type - II distribution. 
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The main objective of this paper is to make comparison between the Bayesian and E-Bayesian method for estimating unknown 
shape parameter of Exponentiated Inverted Weibull Distribution. The Bayes and E-Bayes estimators are obtained by using three 
different loss functions under conjugate prior distribution. Properties of E-Bayes estimates are also studied. 


In this chapter, Bayes and E-Bayes approaches have been used for obtaining the estimates of the unknown parameter of ETWD. In 
section 2, Bayes estimator are derived by using informative prior under Degroot, Al-Bayyati and Minimum Expected loss function. 
In Section 3, Expected — Bayes estimators are derived based on conjugate prior for the parameter of interest. Properties of E- 
Bayesian estimates are discussed in Section 4. 


2. Bayesian Estimation 


In this section, Bayesian estimators of the unknown shape parameter of EIWD are obtained by using gamma informative prior under 
three different loss functions. Let a random sample of size n be drawn from EIWD with pdf given by eq. (1.1), then the likelihood 
function is given by:- 


L(x ,8,8)= [Ti f (9,8) ; 
where x = (X1,X2,X3 Xq4«......X,) and 0 is unknown and £ is known parameter. 
= (6B) [T, xB tY e OZR? (2.1) 


Here, we are using the gamma distribution as an conjugate prior distribution of 6 .The gamma prior distribution with shape 
parameter and scale parameter c and r respectively, has the following pdf 


e(G/e,t) = 91 e-®" 56>0,c,1>0 (2.2) 


This prior distribution was first used by Papadopoulos (1978). 


On combining (2.1) and (2.2), and using Bayes theorem, the posterior density of @ givenx is given by 


n 
n = re 

P(9/x) a ep] | x, 8B FD e- 8 2Xi=1 i Ps re gc-1 eo 

t=1 

= k grte-t 6 [Rix P+ 1] 
where k is independent of 0. 
Thus, Posterior density is given by 
0 (xxi 8 + rn n+c-1 -@ (>? x, 7B + r| 

P(Y/y) = eas 8 e~8 i=1¥i Gere 0 (2.3) 


Tntec 
Bayesian estimate of ® under De-groot loss function using gamma prior 
If 6 is anestimator of @ then the De-groot loss function is given by 


A 6-8 
L(6,6)= =. 
For more details about De-groot loss function, one may refer Degroot (1970), optimal statistical decision. 


The Bayesian estimate using De-groot loss function is given by 





Now, E (275) ree 
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E (7°/x) _ (ntc+1) (ntc) 


(SR x78 + 1)? | 


Bayesian estimator under De-groot loss function is given by 


9. = EM) _ 
Oe E(x) 





(n +c +1) 


acre (2.4) 


Bayesian estimate of 8 under Al-Bayyati new loss function using gamma prior :- 


Al-Bayyati (2002) introduced a new loss function given by 
L(0,6)= 0% (6-0) ;cqeR. 


The Al-Bayyati loss function introduces the additional parameter c, , which assists in determining a flatter loss function and it 
specifically generalizes the Square Error loss function. c, can also be considered the order of weighting of the quadratic 
component. 


E (L(0,6)) = [0 (6-0) * P(9/,) do 


lr (nt+ce+cz) T (n+ c+c,+2) a T (n+ c+c,+1) 


= 62 « 
l(n+c) (xP, xB +r)? Pr (n+ c) (DL x78 + Ya ce (nt c) (DL, x78 + r) 





Ci +1 
Bayes estimator under Al-Bayyati loss function is 


n= & Co ere (2.5) 


RP +H) 
Bayesian estimate of 86 under Minimum Expected loss function using gamma prior :- 
Tummala and Sathe (1978) proposed minimum expected loss function (MELF) which is given as 


(@-6)° Be i : 
ee where @ is anestimator of @. 





L (6,6) = 
a f O=8) . as 
E (L(6,6)) = [—.: P(9/,) de 
0 
_ CORP +) 8? 28 Rix P+ 


(n+c-1) (n+c-2) (n+c-1) 


Bayes estimator under Minimum expected loss function is given by 


~ Be (n+c-2) 
Ome => (2.6) 


n = 
i=1 Xi B+ r) 


3. E-Bayesian Estimates 
According to Han (1997), the prior parameters c andr should be selected to guarantee that the prior g( @ /c, r) in (2.2) be a decreasing 


function of 0. The derivative g( 0 /c, r) with respect to 0 is 


< g(0/c,r) = 7 ge? e- (¢-1) —10] 
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Thus, for 0 < c < 1 andr > 0, the prior distribution g( 8 /c,r) is a decreasing function of 0. 


The E — Bayesian estimate of @ i.e. expectation of the Bayesian estimate of 6 is given by 


ben = E(%/x) 
= Ij Ope 1(O,c,r) * drdc ; 
where Q is the domain of c and r for which the prior density is decreasing in @. Og, is the Bayesian estimate of @ under the 


different loss function. For more details, one can refer Han (2009). 


In order to obtain E-Bayesian estimates of 8, we have to choose prior distribution of hyper parameter c and r. These distributions 
are used to study the impact of different prior distributions on E-Bayesian estimation of @. The following distributions of c and r 
are given by 


m(0,¢,7) = *S :0<c<1,0<r< s (3.1) 
m2 (9,¢,r) = =; 0<c<1,0<r< 5 (3.2) 
m3(8,¢,7) = =; 0<c<1,0<r<s (3.3) 


3(a) E-Bayesian Estimation of 9® under De-groot loss function 


E-Bayesian estimates of @ relative to De-groot loss function based on 11,(0,c,7) , is denoted by Oggn¢ , and is obtained by using 


(2.4) and (3.1). 


(n +c +1) 2(s—r) 
OrepG, = [ E Or x8 + O=x +n ss dr dc 








(s-r) 
ral (n +C +1) (1° (Sh x P+ 1) dr) dc (3.4) 
s___(s-r) aie Xi 
Saactay oF = [loa (SEE-)] (s+ Thaw *) - s (3.5) 
On using (3.5) in (3.4), we get 
ha) ¢ vane i 
9rBpG,= — > 43 4 fliog (eee 1“ )| (s+>f1x; 8)-s } (3.6) 
E-Bayesian estimates of @ relative to De-groot loss function based on 7,(0,c,1r) , is denoted by bene, and is obtained by using 
(2.4) and (3.2). 
A ee Lae 
= * r dc 
mp Jo Jo ORF + Ds 
(6+ 2ie xi 2n+3 
: joa ( x8 )r( re 3.7) 


E-Bayesian estimates of @ relative to De-groot loss function based on 73(0,c,1r) , is denoted by 6repc, and is obtained by using 


(2.4) and (3.3). 
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. (n +c +1) ae d 
EBDG3 = (io go 


12 a— +c bait r 
Ss ( Sos dr) dc(3.8) 


-B 
=a oS SH Ee [log (Hes a )| (3.9) 


(SPL. xi7P + vr) i=1 Xi 


On using (3.9) in (3.8), we get 


ie -B 
Orepa; = f- oa {s = eh liog (Se =e di} dc 


2 
Ss i=1 Xi 


s2 


n xB 
_ (nts) {s ~ yn x8 [tog (GES-}]} (3.10) 


3(b) E-Bayesian estimation under Al-Bayyati loss function. 


E-Bayesian estimates of @ relative to Al-Bayyati loss function based on 71,(9,c,1r) ,is denoted by Ogpap , and is obtained by using 


(2.5) and (3.1). 


(n+ c+cyz) “i 2 (s-r) 


~ 1 
OrpaB, = i iP (xP, xi 78 + r) 52 dr dc 





= {[tog (SEBS) (5 + oe) - 5} «Sn +1 +204) (3.11) 


E-Bayesian estimates of @ relative to Al-Bayyati loss function based on 71,(6,c,r) ,is denoted by Org 4p, and is obtained by using 


(2.5) and (3.2). 


BS _ pl ps (+c+er) 1 
Orpap, = J J. (xx 8 + 1) 2 s dr de 


-B 
= | 108( st Xie 1%i )I* (2n+1+2c, ) (3.12) 


4 xi 8 2s 


E-Bayesian estimates of @ relative to Al-Bayyati loss function based on 713(6,c,r) , is denoted by Ogg AB, and is obtained by using 


(2.5) and (3.3). 


(n+ c+c4) 


a 1 2r 
OrpaBy = fi So Gi ap-Fe 7) * @ dr dc 





2 1 srtyn xj - yh x;,78 
= 24h) (nt coy) (ff Bee Hes ar de (3.13) 
Ec r+ Yh yxiF hx F d _ yn ] St Xie 1%i aa 3 14 
0 (Shix-P +r) r=s—YiLi x? | log m xP (3.14) 


On using (3.14) in (3.13), we have 


as mn .—-B 
Bevan, = SANS — yn xP | log (SHEL) (3.15) 


i=1*i 


s2 
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3(c) E-Bayesian estimation under Minimum Expected loss function. 


E - Bayesian estimates of @ relative to Minimum expected loss function based on 11, (8,c,71), is denoted by Orpy Ff, and is obtained 


by using (2.6) and (3.1). 


a _ pics (n+c-2) 2 (s-r) 
Ozeme, = So So aap dr dc 





=f 2 (n+c—2) Ce (s-r) dr) de 


0 s2 0 (XL, x78 + r) 





n_y,-B a 
= {[tog ( ET ) (s + ver) —s } x = x a 2 (3.16) 


E-Bayesian estimates of @ relative to Minimum expected loss function based on 712(@,c,r) , is denoted by Orpy f, and is obtained 


by using (2.6) and (3.2). 


a 1 fs (n+c-2) 1 
Oreme, = Sy J, (sz i a dr de 


i=1%i * + r+) 


Nop = 
= [log (SPE) | ad (3.17) 


yay 8 2s 
E-Bayesian estimates of @ relative to Minimum expected loss function based on 113(6,c,r) , is denoted by Ore fF, and is obtained 
by using (2.6) and (3.3). 


a 1 fs (n+c-2) 2r 
Greme,= Jy J, Gate ee ae 


i=1 s 


acing = s r+ Cixi? ~Dhaxi? 
=3*J, m+e-2) (J; Gt acFs x) dr )dc 


From equation (3.14), we have 





Q (Tha xiF + 1) WE iF 


Therefore, we have 


- mM .-B 
Brome, = 2{s— Te, xr* | log(*tPz* )Ih G.18) 


s? yas xB 
4. Properties of E-Bayesian estimators . 
In this section, we discuss the relations among different E-Bayesian estimators obtained 


under De-groot loss function , Al-Bayyati loss function and Minimum expected loss function Le. 
Gren; Frpas, 20d Gzemue, (i = 1,2,3) 


Theorem 4.1:- E-Bayes estimator obtained under De-groot loss function follow the following results:- 


1 Orepc, <FEppG, < Feepa,: 
ll. “ lim, OrepG, = lim, OrepG, = “ lim, OrepG;: 
dizi Xi? 7 Lisi %i ? > © ISL ME 08 


Proof:- From eq.(3.6) and (3.7) we get 
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a x n y,-B n y,-B 
Orepc, — FxBdG.= a {log (AES) * ‘aa + 7] = 1} (4.1) 








yy, xj 78 Ss 2 


From equation (3.7) and (3.10), we get 


es ims n_ y.-B n_y.-B 
OrepG,~ F:BpG, = anes {log (ESS * (= + :) = 1} (4.2) 


ny. 
s 4X 





From equation (4.1) and (4.2), we get 


Prep, ~ Peep, = PzBnG,- FEBDG; 


n x.-B n_ x.B 
= 2n+3 {log (Sa ) * (a4 + *) 1} (4.3) 


s rp x7 F s 2 








For -1 < y < 1, we have 


2 3 
In(1 ayaa 











Assuming y = ———~ ; when 0 < —~—- <1 t 
ssuming y Tx e ? when axe , we ge 
ant3 () (s+ ¥7L, x78 ? yr, rhea | 
s og rx 8 s 2 
enh {( = s? + s° (HBS +) -1} 
s Us PRP)” a(R x8)” : ' 
2 3 
mali “| a a a ee at >0 (4.4) 
: 12(DR x8) 6 (DL xi) 


=> Orenc, < Frepc, < Feen«, - 


Proof (ii) :- From equation (4.3), we get 


m 5 ss = 2n+3 (s + yx 8 pe ge 
Orpna, = Prepac, = Prep, = Orppas = 5 {log (ABA Pe | 





n = 
ha iF s 2 
_ im Seppc,~ FEBDG, = <n lim, SrepG,~ Fzevc, = 9 
din1 %i? > © Dini %iP > & 
That is, 7 lim, OrepG, er lim, Oren, Sian lim, OrBDG: 
in Xi F > Dis xP > Zier * PP 7 


Thus, the proof of the theorem is complete. 


Theorem 4.2:- E-Bayes estimator obtained under Al-Bayyati loss function follow the following results:- 


1 Orpas, < Frpas, < Fxpap,: 
ii. lim O nai = lim Onvap = lim as . 
Lh x8 > w * ER xP + 2 Yh x78 + 0 3 


Proof (i) :- From (3.11) and (3.12) 


Grpas, — FrBaB, 





Research Paper 


E-ISSN No : 2454-9916 | Volume: 3 | Issue: 5 | May 2017 





n ,.-B 
= ee flog(1 Sis s ) a (252 + :) -_ 1} 


== 
nx s 2 


From equation, (3.12) and (3.15) we have 


i=1“i 


From (4.5) and (4.6), we get 


Grpas, — Sspas.=FEeBaB, — Fxpas, 


(4.5) 


ot vat n= y.-B 
Orpas, — Fnpap, = Crete) flog (1 + =) - (ee" + :) = 1} (4.6) 




















— (2nt142e1) _s__), (Zax? 1) _ 
= : flog(1 + so) * ( see *) 1} (4.7) 
2 3 
| (1 +55) - (= =S = ) 
OBE TSE 8) SEP amy A aOR 
~ A (2n+1+2c, ) s? s 
ee ' Joo 
EBAB, — VEBAB, ; a (oh x8) o(eh x8) 
> Oren, < Oren, < 6epan,- 
Proof (ii) :- From equation (4.7), we have 
Orpap, = Orpap,= Orpap, = Orpap; 
2n+1+2c s Be Po feck 
5 CAT den) 1) log( 1+—=5 * it Xi + 1 
s aa P s 2 
yn eds OrpaB, = OrpaB, = a gh ops OrpaB, = OrpaB, = 0 
Thus the proof is complete. 
Theorem 4.3:- E-Bayes estimator obtained under Minimum-expected loss function follow the following results:- 
i. 6eeme, < 6epme, < 6reme,- 
ii. lm 6, =. lim: "6 =. lim, 6 
SR ap +00 EBME, et eee EBME> SR xP +00 EBME3 
Proof (i) :- From (3.16) and (3.17) 
6rpme, = 6reme, 
_ 2n-3 s s+2 DR yxi 8) 
aa {log(1 + wat <3] * (ae ) 1} (4.8) 
From equation, (3.17) and (3.18) we have 
=~ =~ 2n-3 s+2 YL x78 


From (4.8) and (4.9), we get 


Oreme, — Seeme,= Feeme, — FeBmes 
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= n ,.-B 
= fog(+ gpg) “(EE = 410 








Ss i=1 Xj 


s s2 s3 


lo 0+ 3)- (5 - SS Se") 
A St ee Ly Pe 


a A 2n-3 s? 33 
OrBABy = OepaB,= a at >0 


i aS = 
s te (ZR ixi78) 6 (DL yx-8)” 
=> Orpas, < Fepap, < Frpas,: 


Proof (ii) :- From equation (4.10), we have 
Oreme, = 6eeme,= 6reme, = 6reme, 


i lim, Oreme, — Srpme, = ae lim, Sepme, — S:eme, = 0 
Diet Xi 7 0 Die Xi sds) 


Thus the proof is complete. 


The part (i) of theorem 4.1, 4.2 and 4.3 shows that with different priors (3.1) — (3.3) of the hyper parameters c and r, the corresponding 
E-Bayes estimates SEBpG, ; 6epan, and Oepme; G=1, 2, 3) are different. The part (ii) of both the theorems shows that SEBDG, ; 


6epan; and Oepme; ; (i=1, 2, 3) are asymptotically equivalent to each other as pare cma — oo, that means SEBDG, (i = 1,2,3) are 
all close to each other when Y_, x;~* is sufficiently largeand Depp; 5 Orem; ;(1=1, 2, 3) are also close to each other . 
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